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ABSTRACT 

We present a computational study of the dynamic behavior of a Ziff-Gulari-Barshad 
model of CO oxidation with CO desorption on a catalytic surface. Our results pro- 
vide further evidence that below a critical desorption rate the model exhibits a non- 
equilibrium, first-order phase transition between low and high CO coverage phases. 
Our kinetic Monte Carlo simulations indicate that the transition process between 
these phases follows a decay mechanism very similar to the one described by the 
classic Kolmogorov-Johnson-Mehl-Avrami theory of phase transformation by nucle- 
ation and growth. We measure the lifetimes of the metastable phases on each side 
of the transition line and find that they are strongly dependent on the direction of 
the transformation, i.e., from low to high coverage or vice versa. Inspired by this 
asymmetry, we introduce a square-wave periodic external forcing, whose two param- 
eters can be tuned to enhance the catalytic activity. At CO desorption rates below 
the critical value, we find that this far-from-equilibrium system undergoes a dynamic 
phase transition between a CO2 productive phase and a nonproductive one. In the 
space of the parameters of the periodic external forcing, this nonequilibrium phase 
transition defines a line of critical points. The maximum enhancement rate for the 
CO2 production rate occurs near this critical line. 



INTRODUCTION 

The study of phase transitions and critical phenomena in nonequilibrium statisti- 
cal systems have recently attracted a great deal of attention due its applications 
in many branches of physics, chemistry, biology, economics, and even sociology [1]. 
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Within this field, surface reaction models have become an archetype for studying 
out-of-equilibrium critical phenomena, and they have been intensely analyzed with 
the purpose of designing more efficient catalytic processes [2]. The Ziff, Gulari, Bar- 
shad (ZGB) model [3] with desorption (the ZGB-k model) [4,5,6] describes kinetic 
aspects of the gas-phase reaction CO+O — > CO2 on a catalytic surface in terms of 
two parameters: the relative partial pressure of CO, y, that represents the proba- 
bility that the next molecule arriving to the surface is a CO, and the desorption 
rate, k, which is related to the probability that an adsorbed CO molecule is desorbed 
without being oxidized. The overall reaction is assumed to occur according to the 
Langmuir-Hinshelwood mechanism, 

CO(g) + E -> CO(a) 

2 (g) + 2E -> 20(a) (1) 
CO(a) + 0(a) -> C0 2 (g) + 2E, 

where E is an empty site on the surface, and (g) and (a) refer to the gas and adsorbed 
phase, respectively. 

SIMULATIONS AND RESULTS 

The ZGB-k model is simulated on a square lattice of linear size L that represents the 
catalytic surface. The kinetic Monte Carlo simulation generates a sequence of trials: 
CO or O2 adsorption with probability 1 — k, and CO desorption with probability 
k. For the adsorption, a CO or 2 molecule is selected with probability y or 1 — y, 
respectively. We calculate the coverages 6*00 and 9o, defined as the fraction of surface 
sites occupied by CO and O, respectively, and the rate of production of CO2, Rco 2 - 
In Fig. [T] we show P(9qo), the probability distribution for 9co vs y, where it is clearly 
seen that at a particular value of y, y2{k), a low 9qo and a high 9q phase coexist. 



P(0 CO ) 




Figure 1: Order-parameter probability distribution, P(9qo), for k = 0.02 and L = 
100. The distribution for the value of y closest to the coexistence value, y2, is shown 
with a bold line. 
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Figure 2: Decay times as functions of y when the system evolves toward the low CO 
coverage region (r d ), and when it evolves toward the high coverage region (r d ), shown 
for k = 0.02 and L = 100. Here and elsewhere in this paper, time is measured in 
Monte Carlo steps per site (MCSS). 

A finite-size-scaling analysis of the statistical fluctuations of the CO coverage gives 
strong evidence that below the critical value of k the model exhibits a first-order, 
nonequilibrium phase transition between low and high CO coverage phases with the 
same characteristics as a first-order equilibrium phase transition [4,5,7]. 

We also measured the metastable lifetimes associated with the transition from 
the low CO coverage phase to the high CO coverage phase and vice versa. The 
system-size dependence of the decay times strongly suggests that the system follows 
a decay mechanism very similar to the one described by the classic Kolmogorov- 
Johnson-Mehl-Avrami theory of phase transformation by nucleation and growth near 
a first-order equilibrium phase transition, with well-defined single-droplet and mul- 
tidroplet regimes [8,9]. In this system, the desorption parameter and the distance to 
the coexistence curve play the roles of the temperature and the supersaturation or 
overpotential, respectively. Near the coexistence curve the decay times are inversely 
proportional to l/L, and the decay mechanism consists of the nucleation and growth 
of a single droplet of the stable phase. Far from the coexistence curve, the decay times 
are independent of the system size, and the decay proceeds by random nucleation of 
many droplets of the stable phase [5]. 

We found that the lifetimes strongly depend on k and on the direction of the 
process; the mean decontamination time (r^) (from high to low CO coverage) being 
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Figure 3: (a) Applied periodic pressure of CO, y(t), that takes the values y\ = 0.5 and 
yh = 0.5346 during the time intervals = 10 and t p = 50, respectively, (b) Response 
of the production rate to the applied pressure given in (a) for L = 100 with k = 0.02. 
The dotted line marked (r) indicates the long-time average of the period-averaged 
CO2 production rate r, while the dotted line marked (-Rco 2 )max marks the maximum 
average C0 2 production rate for constant y = yz(k). 



different from the mean poisoning time (r p ) (from high to low CO coverage). At 
comparable distances from the coexistence curve, (rA (r p ), as seen Fig- El Since 
several experiments indicate that it is possible to increase the efficiency of a catalytic 
process by subjecting the system to periodic forcing [10], we decided to exploit the 
asymmetry between the decay times by subjecting the system to a periodic variation 
of the external pressure with periods related to the decay times in each direction. We 
therefore select a square-wave periodic variation of the CO pressure, y(t), that in a 
period T = td + t p takes the values, yi (located below the transition pressure) during 
the time interval td and yt (located above the transition pressure) during the time 
interval t p , as can be seen in Fig. |3fa). 

We found that the times that the driving force spends in the low and high coverage 
regions, t d and t p , respectively, can be tuned for each set of yi and yh to increase 
the productivity of the system. In Fig. Efb) it is seen that the CO2 production 
rate exhibits an oscillatory behavior in response to the periodic pressure shown in 
Fig. Eta)- The period-averaged value of the C0 2 production rate Rco 2 i defined as 

r = ± f Rco a (t)dt, (2) 
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Figure 4: Long-time average of the period-averaged rate of C0 2 production, (r), 
shown vs t d for two values of t p and L = 100; (a) with yi = 0.52, y h = 0.535, and 
k = 0.01, and (b) with yi = 0.52, y^ = 0.553, and k = 0.04. Only for the lower value 
of k does the system clearly present two dynamic phases: one with (r) w and the 
other with (r) > 0. 



plays the role of the dynamic order parameter [6]. For the parameters selected in 
Fig. El the long-time average of r, (r) takes a value that is about 7% higher than the 
maximum average CO2 production rate for constant y. It is likely that more careful 
tuning of the parameters could further improve the degree of enhancement. 

We also found that, for sufficiently low values of k, the driven system undergoes 
a dynamic phase transition between a dynamic phase of high CO2 production, r > 0, 
and a nonproductive one, r w 0, as can be seen in Fig. 0] The distinction between 
these phases disappears for desorption rates above the critical value [6]. 

A detailed finite-size scaling analysis indicates that for small values of k, the 
measure of the fluctuations of the order parameter, 

X L = L 2 [(r 2 ) - (r) 2 ] , (3) 

diverges as a power law with the system size, X™ ax w L 7 ^ with exponent 7/1/ = 
1.77 ± 0.02, while moments of the order parameter at the transition point decay as 
{r n )i m L~ n ^l v with (3/v = 0.12 ±0.04. These exponent ratios, together with general 
symmetry arguments, give reasonable evidence that this far-from equilibrium phase 
transition belongs to the same universality class as the equilibrium Ising model [6]. 
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Figure 5: Long-time average (r) of the period-averaged CO2 production rate r, shown 
as a density plot vs t p and for y\ = 0.51, = 0.535, k = 0.01, and L = 100. 

The long-time average production rate (r) is shown in a density plot vs t p and t& in 
Fig. |SJ The line of critical points appears as the sharp boundary of the black, low- 
production region in the lower right-hand part of the figure. The region of maximum 
average production is seen to lie very close to the critical line on its high-production 
side. We believe this observation offers a clue to understanding the physical reason 
for the enhancement. Most likely, the long-range critical correlations associated with 
the critical cluster that develops for parameter values near the critical line produce 
a high density of sites where CO molecules would be situated next to O atoms, in 
positions conducive to rapid oxidation and desorption of the produced C0 2 . 

CONCLUSION 

In this paper we have summarized results of large-scale kinetic Monte Carlo simula- 
tions of the ZGB-k model of heterogeneously catalyzed CO oxidation with CO desorp- 
tion. The results of the simulations were analyzed using finite-size scaling methods. 
We found that there are strong similarities between the dynamics of metastable decay 
in this far-from equilibrium, non-Hamiltonian system and the well-known behavior 
of Hamiltonian systems. From a theoretical point of view, these similarities could 
lead to significant advancement in our understanding of the dynamics of far-from- 
equilibrium systems. From a practical point of view, our results can be exploited to 
develop novel ways to increase the efficiency of catalytic reactions. 
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